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1.  BACKGROUND 

Gas  guns  are  used  at  the  Harry  Diamond 
Laboratories  (HDL)  to  accelerate  ordnance 
items  to  velocities  typical  of  mortar,  recoilless 
rifles,  and  artillery  projectiles,  that  is,  400  to 
4000  ft/s  (120  to  1200  m/s).  Test  considera¬ 
tions  often  dictate  achievement  of  the  required 
velocity  with  the  lowest  acceleration  (g)  force 
exerted  on  the  projectile  while  in  the  gun.  This 
dictum  necessitates  the  use  of  low  driving 
pressures,  and  therefore  long  guns  are  needed 
to  achieve  the  required  velocities.  These  guns 
are  installed  in  special  purpose  laboratory 
space,  and  the  maximum  gun  length  is  limited 
by  the  dimensions  of  the  building.  The 
characteristics  of  the  guns  and  the  governing 
theories  must  be  known  to  configure  the  guns 
so  as  to  maximize  the  projectile  velocity  for  the 
applied  acceleration. 

The  gun  used  most  often  at  HDL  has  a 
4-in.  bore  and  is  the  subject  of  this  experimen¬ 
tal  and  theoretical  investigation.  A  large  por¬ 
tion  of  the  gun  is  located  in  a  5-ft  (1.5-m)-high 
tunnel  that  runs  the  width  of  the  building  that  it 
occupies.  Several  meters  of  each  end  of  the 
gun  extend  into  a  "breech”  room  and  a  "muz¬ 
zle"  room.  This  gun  can  be  configured  in  two 
ways  (fig  1).  For  many  tests,  velocities  of  400 
to  800  ft/s  (120  to  240  m/s)  are  adequate. 
These  are  readily  achievable  in  a  1-atm  con¬ 
figuration.  In  the  1-atm  gun,  a  pin  restrained 
projectile  seals  one  end  of  a  barrel,  a 
diaphragm  seals  the  other  end,  and  the  barrel 
is  evacuated  (fig.  la).  Removal  of  the  pin  per¬ 


(a) 


Figure  1.  Four-inch  gun:  (a)  one-atmosphere  and  (b) 
pressurized  configurations. 


mits  atmospheric  pressure  to  accelerate  the 
projectile  down  the  barrel  and  through  the 
diaphragm  When  tests  require  higher 
velocities,  a  pressurized  driver  is  required.  The 
gun  contains  provision  for  a  second  pin 
downstream  from  the  1-atm  pin.  The  breech 
end  of  the  gun  is  sealed,  and  the  section  be¬ 
tween  this  end  of  the  gun  and  the  projectile 
can  be  pressurized  (fig.  1b).  This  configuration 
uses  part  of  the  gun  length  as  driver  and 
therefore  shortens  the  length  of  the  gun 
available  as  barrel. 

Consider  the  four  pressure-driven  gun  con¬ 
figurations  shown  in  figure  2.  In  the  uppermost 
gun  (fig.  2a),  the  driver  gas  has  no  volume,  and 
therefore  the  projectile  can  achieve  no  veloci¬ 
ty.  The  projectile  moves  with  increasing  speed 
as  the  pin  is  located  farther  downstream  (fig. 
2b).  Further  movement  of  the  pin  (fig.  2c) 
shortens  the  barrel  length  (the  distance  over 
which  the  projectile  accelerates)  and  thus  off¬ 
sets  the  gain  in  chamber  volume.  Finally,  as 
the  pin  location  nears  the  muzzle  end  of  the 
gun  (fig.  2d),  the  projectile  has  little  distance 
over  which  to  accelerate,  and  the  muzzle 
velocity  approaches  zero. 

The  problem  is  to  locate  the  pin  at  the  posi¬ 
tion  that,  for  any  specified  velocity,  requires 
the  lowest  gas  pressure.  This  is  equivalently 
stated  as  the  pin  position  that  maximizes  the 
projectile  muzzle  velocity  for  a  given  gas 
pressure.  There  are  a  limited  number  of  pin 
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Figure  2.  Four  locations  for  release  pin  of  constant 
length  gun:  (a)  at  breech,  (b)  near  breech,  (c)  near 
muzzle,  and  (d)  at  muzzle. 


positions  (breech  locations)  in  the  4-in  gun 
because  of  the  way  that  the  gun  is  con¬ 
structed.  This  report  shows  that  no  one  pin 
position  is  optimal  for  all  velocities  and, 
because  it  is  impractical  to  relocate  the  breech 
for  different  tests  the  solution  to  the  problem 
becomes  somewhat  subjective. 


2.  CONSTRUCTION 

The  gun  used  in  this  investigation  is  con¬ 
structed  from  eight  12-ft  (4-m)  lengths  of 
90-mm  gun  barrels.  The  barrels  have  been 
bored  to  a  4-in.  (10-cm)  inside  diameter  (i.d.), 
turned  down  to  a  6-in.  (1 5-cm)  outside  diameter 
(o.d.),  and  flanged.  The  barrel  sections  are 
always  evacuated,  even  when  a  pressure 
driver  is  used.  The  joints  are  sealed  with  cop¬ 
per  gaskets.  The  muzzle  section  of  the  gun  (fig. 
3)  attaches  to  the  last  barrel  and  contains  the 
fittings  for  the  vacuum  connections.  This  sec¬ 
tion  has  a  large  internal  toroidal  reservoir  to 
accept  any  gas  that  is  driven  ahead  of  the  pro¬ 
jectile  and  thereby  to  minimize  the  pressure 
rise  ahead  of  the  projectile.  This  section  con¬ 
tains  also  the  muzzle  seal,  which  is  a  0.002-in 


Figure  3.  Muzzle  section  of  4-in.  gun. 


(50-^m)  Mylar  sheet  clamped  in  a  door  by  two 
screw  latches.  The  Mylar  is  located  on  a  roller 
above  the  door,  and  a  new  section  of  the  sheet 
is  pulled  down  and  reclamped  after  each  shot. 

The  release  pin  assembly  used  in  the 
1-atm  configuration  (fig  4)  attaches  to  the  en¬ 
trance  of  the  first  barrel.  It  is  a  short  gun  exten¬ 
sion  with  an  internal  circumferential  0-ring 
groove  to  seal  to  the  girth  of  the  projectile  and 
a  pneumatically  operated,  electrically 
controlled  pin  to  restrain  the  projectile.  Also 
visible  is  a  bac1  jiate,  which  acts  as  a  backstop 
to  capture  the  projectile  if  it  is  propelled  out  of 
the  rear  of  the  gun  in  the  event  that  the  muzzle 
diaphragm  is  accidentally  ruptured.  To  the  left 
of  the  release  pin  assembly  is  a  side  access 
port,  which  remains  covered  in  the  1-atm  con¬ 
figuration  and  is  used  only  in  the  pressurized 
configuration. 

In  the  pressurized  configuration,  the 
breech  end  of  the  gun  is  sealed  by  a  closure 
section  located  in  the  side  access  port  in  the 
first  barrel  section  (fig.  5)  The  closure  section 
contains  pressurization  and  vacuum  ports. 

The  projectile  restraining  section  used  in 
the  pressurized  configuration  (fig.  6)  is  1  ft  (30 
cm)  long;  has  a  side  access  panel  8  in.  long; 
and  contains  a  pnuematically  operated,  elec¬ 
trically  controlled  pin  and  the  required  0-ring. 
The  projectile  restraining  section  can  be 
placed  between  any  two  sections  of  the  gun 
and  therefore  can  be  placed  at  one  of  seven 
possible  locations  along  the  barrel.  When  the 
1-atm  configuration  is  used,  the  projectile 
restraining  section  0-ring  is  removed,  and  the 
pin  is  withdrawn  to  allow  free  projectile 
passage  through  the  projectile  restraining  sec¬ 
tion 
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Figure  4.  One-atmosphere  configuration  breech  of  4-in.  gun. 


Figure  5.  Closure  section. 
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Figure  6.  Projectile  restraining  section:  (a)  open  and  (b)  closed. 
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The  drawing  in  figure  7  shows  the  4-in.  gun 
with  the  projectile  restraining  section  mounted 
between  the  first  two  barrel  sections.  This 
assembly  provides  easy  access  to  the  projec¬ 
tile  restraining  section  because  it  is  located  in 
the  breech  room.  Currently,  the  projectile 
restraining  section  is  mounted  between  barrel 
sections  3  and  4  and  is  located  in  the  gun  tun¬ 
nel.  This  location  is  somewhat  inconvenient 
because  of  the  low  ceiling 


ItitoNi.  I? 

I  flair  j  a1  =r  n1  ift-t* ) 

Figure  7.  Isometric  side  view  of  4-in.  gun. 

3.  VELOCITY  MEASURING  INSTRUMEN¬ 
TATION 

Muzzle  speed  is  obtained  by  measurement 
of  the  time  of  flight  between  two  “light 
screens"  located  12.00  in.  (30.48  cm)  apart 
near  the  end  of  the  last  gun  barrel  section.  All 
the  experimental  projectile  muzzle  velocity 
data  presented  below  are  computed  by  using 
these  devices.  For  theoretical  purposes,  the 
end  of  the  gun  is  assumed  to  be  halfway  be¬ 
tween  the  two  light  screens. 

Each  light  screen  consists  of  a  light 
source  operated  by  direct  current  (dc),  a 
photodetector,  and  an  amplifier  circuit.  The 
source  and  the  detector  face  each  other  along 
a  diameter  of  the  gun  tube.  When  the  projectile 
interrupts  the  light  falling  on  the  detector,  a 
trigger  signal  is  generated. 

The  light  source  consists  of  an  incandes¬ 
cent  lamp  mounted  behind  an  achromatic  lens, 
which  images  the  lamp  filament  about  200  mm 
in  front  of  the  lens.  A  lens  in  the  detector 
assembly  refocuses  the  light  onto  a  small 
piece  of  frosted  glass  positioned  immediately 
in  front  of  the  photodetector.  With  this  arrange¬ 
ment,  the  photodetector  amplifier  output  jumps 
from  2  to  18  V  when  the  light  beam  is  inter¬ 


rupted.  The  rise  time  of  each  detector  and 
amplifier  combination  is  15  ns,  with  less  than 
1-ms  difference  between  the  two  light  screens. 
Therefore,  time  interval  measurements  are  ac¬ 
curate  to  within  1  /is;  positioning  of  the  light 
screens  is  accurate  to  within  0.01  in.  (0.25  cm); 
hence,  the  accuracy  of  the  velocity 
measurements  is  approximately  0.4  ft/s  (12 
cm/s)  at  400  ft/s  (1 20  m/s)  and  about  6  ft/s  (1 .8 
m/s)  at  2500  ft/s  (750  m/s). 

Previous  setups  with  finer  beams  and 
therefore  finer  distance  resolution  proved  to  be 
sensitive  to  debris  or  vibrations  or  both  and  oc¬ 
casionally  triggered  early.  The  present  ar¬ 
rangement  has  proven  to  be  highly  reliable. 

Streak  photography,  which  is  used  to 
photograph  the  projectile  shortly  after  it  leaves 
the  muzzle,  provides  a  further  check  on  muzzle 
velocity.  The  streak  method  (used  for  muzzle 
velocities  less  than  1500  ft/s  or  300  m/s) 
agrees  with  the  above  method  to  better  than  1 
percent. 


4.  OTHER  MEASUREMENTS 

Other  measurements  that  enter  into  the 
calculation  of  the  initial  acceleration  are  pro¬ 
jectile  weight,  gas  pressure,  gun  length,  and 
gun  diameter.  By  far,  the  least  accurate  of 
these  measurements  is  the  gas  pressure, 
which  is  good  to  only  5  percent;  hence, 
measurements  of  the  initial  acceleration  are 
good  to  about  5  percent. 

Certain  gas  parameters  are  used  in  the 
theory.  The  gas  properties  assumed  are  the 
properties  of  a  pure  gas  at  20  C.  The  room 
(driver)  temperature  was  constant  to  about  3  C, 
which  has  little  effect  on  these  properties.  For 
nitrogen  experiments,  no  special  care  was 
taken  to  assure  driver  gas  purity.  It  was  as¬ 
sumed  that  the  air  initially  in  the  driver  section 
was  sufficiently  close  in  properties  to  nitrogen 
so  as  not  to  influence  the  results.  For  helium 
shots,  a  purging  technique  was  implemented 
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consisting  of  three  series  of  evacuation  of  the 
driver  section  to  less  than  5  psia  (34.5  kPa  ab¬ 
solute)  followed  by  helium  pressurizations  to 
50  psig  (345  kPa  gauge)  (except  for  the  final 
pressurization,  which  was  to  the  required  firing 
pressure). 

5.  THEORY 

The  problem  of  optimization  of  pin  position 
in  the  4-in.  gun,  along  with  a  more  general 
problem  that  includes  diameter  changes  in  the 
driver  section,  has  been  analyzed  by  Seigel.1 
Seigel’s  report  contains  an  analysis  of  the  four 
active  HDL  gas  guns  and  two  future  guns.  The 
aim  of  his  investigation  is  to  configure  the  guns 
so  as  to  maximize  the  projectile  velocity  for  a 
specified  maximum  projectile  acceleration. 
The  aim  of  HDL’s  investigation  is  stated  in  sec¬ 
tion  1. 

5.1  Isentropic  Expansion 

It  is  instructive  to  start  the 
theoretical  discussion  with  a  problem  that  is 
analytically  tractable,  adds  insight  to  the 
results,  and  is  shown  to  be  valid  for  certain  test 
conditions.  Isentropic  expansion  of  a  driver 
gas2  yields  an  approximate  prediction  of  gun 
muzzle  velocity  (Lagrange  approximation) 
when  the  gas  sound  speed  is  much  larger  than 
the  projectile  velocity  (U),  that  is,  small  Mach 
number.  It  is  assumed  that  the  driver  section 
diameter  is  the  same  as  the  barrel  diameter. 

U2  =  K[NZ(1  -Z2,n)/2]  ,  (1) 

where 

K  =  2APL/M  .  (ia) 

Z  =  X/L  .  (1b) 

N  is  the  degrees  of  freedom  of  the  driver  gas, 

1 Arnold  E.  Seigel,  Performance  Calculations  and  Op¬ 
timization  of  Gas  Guns,  Chevy  Chase,  MD,  HDL- 
CR-8 1-723-1  (May  1981). 

2J.  Corner,  Theory  of  the  Interior  Ballistics  of  Guns,  John 
Wiley  &  Sons,  New  York  (1950). 


A  is  the  cross-sectional  area  of  the  gun  barrel, 
P  is  the  initial  gas  pressure  in  the  driver  reser¬ 
voir,  L  is  the  total  gun  length,  M  is  the  projectile 
mass,  Z  is  the  relative  position  of  the  pin 
measured  from  the  breech  end  of  the  barrel, 
and  X  is  the  length  of  the  driver  reservoir. 

The  square  root  of  K  is  the  velocity 
(U totai)  that  would  be  achieved  in  a  constant  ac¬ 
celeration  gun  of  the  total  gun  length. 

U,o,al  =  K1'2.  (ic) 


The  quantity  in  brackets  in  equation  (1)  is 
always  less  than  1  and  therefore  represents  a 
decrease  from  Ut0!a,. 

Since  the  ratio  of  barrel  length  (Y)  to 
total  gun  length  is  1  -  Z, 


U barrel  =  U total  0  —  Z) 1/2  ■  (Id) 

Uba(re,  is  the  velocity  that  would  be  attained  if 
the  acceleration  were  constant  over  the  barrel 
length. 

The  optimum  pin  position  (Zopt)  is 
found  by  differentiating  U  with  respect  to  Z  and 
setting  the  result  to  zero.  These  actions  result 
in 


Zop,  =  (1  +2/NrN,z  , 

(2) 

ft  max  =  ft  lotal  [Zopt  /(I  +  2/N)]1 

(3) 

Umax  =  U barrel  [N/(N  +  2)],/2 

*[zwi(i-zop[)V12  ■ 

(4) 

Table  1  lists  results  from  equations 
(2)  to  (4).  Note  these' 

a  The  greater  the  degrees  of 
freedom  of  the  gas,  the  shorter  the  driver  and 
the  higher  the  expected  velocity. 
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b  When  the  pin  is  located  about  40 
percent  of  the  way  down  the  length  of  the  gun 
barrel,  we  should  expect  about  55  percent  of 
the  velocity  achievable  if  the  entire  gun  length 
were  available  for  constant  acceleration 

c.  Velocity  based  on  the  actual  bar¬ 
rel  length  indicates  that  the  projectile  can 
never  achieve  better  than  three-quarters  of  the 
constant  acceleration  velocity. 

TABLE  1.  FEATURES  OF  OPTIMIZED  LAGRANGE 
GUN 


N 

^opt 

^max^total 

Umax^U  barrel 

3 

0.465 

0.528 

0.722 

5 

0.431 

0.555 

0.736 

Infinite 

0.368 

0.606 

0.763 

5.2  Computer  Analyses 

Theories  that  account  for  inertial  ef¬ 
fects  of  the  gas  were  programmed  in  VAX-11 
BASIC  and  are  treated  in  the  following  sec¬ 
tions.  The  programs  are  included  in  appendix  A 
and  are  referred  to  as  used  in  this  paper.  For 
ease  of  presentation,  the  results  are  in 


graphical  form.  The  results  obtained,  although 
sometimes  presented  as  dimensionless 
numbers,  are  not  generally  applicable.  Guns  of 
different  lengths  yield  different  curves  as  func¬ 
tions  of  the  same  dimensionless  numbers. 
Such  curves  can  be  obtained  only  by  changing 
the  initial  conditions  in  the  programs. 

5.3  Pidduck-Kent  Special  Solution 

An  important  analytic  treatment  is 
the  Pidduck-Kent  solution,  which  assumes  that 
through  a  number  of  reflections  between  the 
rear  of  the  driver  and  the  projectile,  the  flow 
has  reached  a  steady  pressure  profile,  which 
decreases  in  amplitude  with  time.  The  equa¬ 
tions  governing  this  solution  are  included  in  the 
program  PDQ  and  are  not  repeated  here. 

In  the  following  graphs,  projectile 
velocity  is  presented  as  a  function  of  relative 
pin  position,  that  is,  driver  section  length  per 
total  gun  length.  The  parameter  is  the  initial  ac¬ 
celeration,  which  is  the  same  as  the  maximum 
acceleration.  Figure  8  presents  results  for  a 
helium  driver.  The  initial  acceleration  is  200  g 
for  the  lowest  curve  and  3000  g  for  the  highest 
curve,  and  intermediate  curves  are  in  200-g  in¬ 
crements.  Figure  9  shows  the  same  results  for 
nitrogen 
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Figure  8.  Pidduck*Kent  theory,  helium,  200  to  3000  g. 


Figure  9.  Pidduck-Kent  theory,  nitrogen,  200  to  3000  g. 
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Note  that  the  relative  pin  position  for 
the  velocity  peaks  for  the  lowest  curves  in 
each  figure,  0.467  and  0.417,  fall  near  those 
predicted  by  the  Lagrange  theory  in  table  1. 
However,  as  the  initial  acceleration  increases, 
the  peaks  shift  to  smaller  values  of  relative  pin 
position.  The  better  agreement  for  the  helium 
data  is  attributed  to  the  assumption  that  a 
weightless  gas,  inherent  in  the  Lagrange 
theory  (which  postulates  infinite  sound  speed), 
is  better  satisfied;  the  Mach  number  for  the 
200-g  helium  case  is  only  0.15,  but  for  the 
nitrogen  case  it  is  up  to  0.5.  The  lowest  curves 
(200  g)  also  predict  velocities  of  576  and  582 
ft/s  (173  and  175  m/s)  for  helium  and  nitrogen, 


that  is,  higher  for  nitrogen!  However,  as  the  in¬ 
itial  acceleration  increases,  the  higher  sound 
speed  of  helium  becomes  of  significance,  and 
higher  projectile  velocities  are  obtained  with 
the  monatomic  gas. 

Figures  10  and  11  present  data  for 
initial  accelerations  of  3000  to  15,000  g  in 
steps  of  2000  g.  Although  accelerations  of 
more  than  3000  g  are  rarely  anticipated  for  the 
particular  application  of  this  investigation,3 
these  results  are  included  for  completeness. 

3Herbert  D  Curchack,  An  Artillery  Simulator  for  Fuze 
Evaluation,  Harry  Diamond  Laboratories  HDL-TR-1330 
(November  1966). 
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Figure  10  Pidduck-Kent  theory,  helium,  3000  to  15,000  g. 
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Figure  11  Pidduck-Kent  theory,  nitrogen,  3000  to  15,000  g. 


5.4  Constant  Diameter,  Infinite 
Length  Driver  Theory 

If  the  driver  is  assumed  to  be  the 
same  diameter  as  the  barrel  and  infinite  in 
length,"  then  program  CONDIA  is  applicable. 
Although  the  governing  equation  can  be 
deduced  from  the  program,  it  is  camouflaged 
by  changes  in  notation  and  is  repeated  here. 

PAY/Ma2  =  [u-N/(N  +  1)](1  —  u/N)*N+1>  (5) 

+  N/(N  +  1), 

where  a  is  the  initial  sound  speed  in  the  driver 
reservoir  and  u  is  the  projectile  Mach  number 
based  on  initial  sound  speed.  This  theory  ap¬ 
plies  where  the  driver  “appears”  infinite.  That 

4 Arnold  E.  Seigel,  The  Theory  of  High  Speed  Guns,  Naval 
Ordnance  Laboratory,  White  Oak,  MD,  AGARDograph  91 
(I May  1965),  eg  (12.2). 


is,  when  the  first  expansion  wave  (of  an  expan¬ 
sion  fan)  leaves  the  projectile  at  the  onset  of 
motion,  the  wave  is  reflected  from  the  rear  of 
the  gun,  but  does  not  overtake  the  projectile 
before  the  projectile  exits  from  the  muzzle  The 
resultant  curves  are  plotted  in  figures  12  to  15 
for  helium  and  nitrogen  and  for  the  lower  and 
higher  initial  acceleration  ranges. 

The  diamonds  on  the  curves  repre¬ 
sent  the  projectile  starting  location  for  which 
the  first  reflected  wave  overtakes  the  projectile 
just  as  it  leaves  the  muzzle  They  were 
obtained  from  program  CATCH.  The  assump¬ 
tion  of  infinite  length  driver  is  valid  for  all  points 
to  the  right  of  the  diamonds,  and  the  validity 
decreases  to  the  left. 

These  curves  do  not  have  theoretical 
maxima  as  in  the  Pidduck-Kent  and  Lagrange 
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cases.  The  optimum  pin  location  from  these 
curves  is  at  a  relative  pin  position  of  zero  (  Y  = 
L),  which  is  in  a  region  in  which  the  theory  does 
not  hold.  (The  driver  is  of  zero  length,  not  in¬ 
finite  length.)  Therefore,  a  determination  of  op¬ 
timum  pin  position  from  these  curves  is  not 
possible.  However,  the  first  waves  that  over¬ 


take  the  projectile  are  weak  expansion  waves, 
and  theoretical  results  that  include  effects  of 
these  waves  would  follow  infinite  driver  results 
for  some  distance  to  the  left  of  the  diamond. 
Based  on  this  consideration,  the  curves  would 
continue  to  rise  for  some  distance  to  the  left  of 
the  diamonds,  and  optimum  relative  pin  posi¬ 
tion  would  be  to  the  left  of  the  diamonds. 
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Figure  12.  Infinite  driver  theory,  helium,  200  to  3000  g. 
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Figure  13.  infinite  driver  theory,  nitrogen  200  to  3000  g. 
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Figure  14.  Infinite  driver  theory,  heiium,  3000  to  15.000  g. 
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Figure  15.  Infinite  driver  theory,  nitrogen,  3000  to  15,000  g. 


5.5  Combined  Pidduck-Kent  and  Infinite 
Driver  Theories 

The  results  presented  in  figures  16  to 
19  were  obtained  from  program  CDPK,  which, 
combines  the  results  of  Pidduck-Kent  for  small 
values  of  relative  pin  position  with  infinite 
driver  results  for  large  values.  These  were  ob¬ 
tained  by  comparing  the  predicted  velocities 
from  the  two  theories  starting  at  a  relative  pin 
position  of  zero  and  using  the  Pidduck-Kent 
result  as  output  so  long  as  it  predicted  a  veloci¬ 
ty  lower  than  the  infinite  driver  results.  Once 
the  infinite  driver  theory  predicted  lower  veloci¬ 
ty,  infinite  driver  results  were  used  for  the  re¬ 
mainder  of  the  run,  even  though  for  higher 
values  of  relative  pin  position  Pidduck-Kent 
again  predicted  lower  muzzle  velocity.  Infinite 
driver  results  were  used  because  infinite  driver 
theory  is  more  appropriate  in  this  region,  The 


point  at  which  the  two  theories  intersect  is 
marked  with  a  plus  sign.  The  diamond 
represents  the  projectile  starting  location  for 
which  the  first  reflected  wave  overtakes  the 
projectile  just  as  it  leaves  the  muzzle 

Generally,  for  low  initial  accelera¬ 
tions,  the  transition  from  one  theory  to  the 
other  is  smooth,  and  the  first  reflection  points 
are  on  or  very  close  to  the  Pidduck-Kent  por¬ 
tion  of  the  curve.  For  higher  accelerations,  the 
transition  becomes  more  abrupt,  and  the  first 
reflection  points  lie  above  the  Pidduck-Kent 
curves.  Since  infinite  driver  theory  (sect.  5.4) 
implies  that  the  maximum  is  to  the  left  of  the 
diamond,  the  conclusion  drawn  from  these 
figures  is  that  Pidduck-Kent  underestimates 
muzzle  velocity  in  the  transition  region,  that  is, 
where  the  first  reflection  overtakes  the  projec¬ 
tile  near  the  muzzle.  Since  Pidduck-Kent 
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assumes  a  gas  pressure  distribution  that  is 
achievable  only  after  several  reflections  and 
infinite  driver  assumes  no  reflections  from  the 
projectile,  this  divergence  is  not  unexpected. 
To  use  these  graphs  to  predict  optimum 
relative  pin  position  now  becomes  difficult, 
especially  for  nitrogen.  Although  the  curves 
have  well-defined  peaks,  these  peaks  are  in 
the  region  where  the  wrong  theory  defines  the 
peaks.  To  reiterate,  the  “real”  peak  must  be 
someplace  to  the  left  of  the  diamonds. 

Another  consideration  in  predicting 
velocity  is  associated  with  the  very  low  Mach 
number  region.  For  example,  consider  low  ini¬ 
tial  accelerations  with  helium.  The  sound 
speed  is  sufficiently  high  so  that  many  wave 
reflections  occur  while  the  projectile  moves 
only  slightly.  This  situation  leads  to  a  nearly 


uniform  gas  pressure  between  the  projectile 
and  the  end  of  the  breech.  When  the  projectile 
remains  at  a  low  Mach  number,  the  pressure 
drops  uniformly  in  the  driver  section  as  the  pro¬ 
jectile  moves  and  the  volume  increases.  Under 
this  condition,  Lagrange  results  should  be 
more  appropriate  than  Pidduck-Kent  results. 


The  net  result  of  these  theoretical 
considerations  is  that  the  velocity  peak  moves 
toward  lower  relative  pin  position  with  increas¬ 
ing  initial  acceleration,  the  peak  is  somewhere 
to  the  left  of  and  near  the  diamonds,  and  the 
peaks  differ  in  amplitude  and  location  for 
helium  and  nitrogen.  However,  experimental 
results  must  be  evaluated  to  verify  the  ap¬ 
propriateness  and  the  range  of  any  of  the 
theories  presented. 


Figure  16.  Pidduck-Kent  and  infinite  driver  theories  combined,  helium,  200  to 
3000  g. 
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Figure  17  Pidduck-Kent  and  infinite  driver  theories  combined,  nitrogen 
200  to  3000  g. 
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Figure  18  Pidduck-Kent  and  infinite  driver  theories  combined,  helium,  3000  to  15,000  g. 
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Figure  19.  Pidduck-Kent  and  infinite  driver  theories  combined,  nitrogen,  3000 
to  15,000  g. 


6.  EXPERIMENTAL  RESULTS 
6.1  Procedure 


Many  of  the  results  included  in  what 
follows  were  made  during  routine  operation  of 
the  gun.  The  projectiles  were  not  specially 
prepared  for  optimum  gun  performance  (fig. 
20).  Several  projectiles  were  used  The  most 
critical  tolerances  for  these  experiments, 
those  associated  with  projectile  diameter  and 
the  distance  between  the  two  bore  riders 
(wheel  base),  were  only  grossly  controlled.  The 
projectile  diameter  ranged  from  0.004  to  0.015 
in.  (102  to  381  jim)  less  than  the  gun  i  d.  The 
wheel  base  was  typically  1  to  1 .5  calibers.  The 
center  of  mass  of  the  projectile  typically  was 
near  or  aft  of  the  rear  bore  rider. 


6  2  Data 

Experiments  have  been  run  with 
nitrogen  with  the  projectile  restraining  section 
at  two  positions  and  with  helium  with  the  pro¬ 
jectile  restraining  section  at  one  of  these  posi¬ 
tions.  Maximum  pressure  was  300  psia  (2070 
kPa  absolute).  Projectile  weights  varied  from 
0.5  to  5  lb  (0.25  to  2.5  kg)  and  covered  the 
velocity  range  from  500  to  2600  ft/s  (150  to  780 
m/s).  The  first  series  of  experiments  had  the 
projectile  restraining  section  at  the  end  of  the 
first  breech  section;  this  was  the  easiest  sec¬ 
tion  to  access  Sixty-nine  shots  were  fired  us¬ 
ing  nitrogen  with  the  projectile  restraining  sec¬ 
tion  located  here.  Subsequently,  229  tests  with 
nitrogen  and  94  tests  with  helium  were  run  with 
the  projectile  restraining  section  between  gun 
sections  3  and  4. 
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Figure  20.  Typical  projectiles  used  in  study. 


Each  set  of  test  data  is  plotted  along 
with  predictions  from  each  of  the  three 
theories  (fig  21  to  23).  The  theoretical  predic¬ 
tions  were  obtained  from  the  program 
THEORY. 

If  the  gun  behaved  as  a  constant  ac¬ 
celeration  gun,  then  the  velocity  would  be  pro¬ 
portional  to  the  square  root  of  the  acceleration 
Therefore,  the  x-axis  in  the  figures  was  chosen 
to  be  the  square  root  of  the  initial  acceleration 
rather  than  the  initial  aceleration.  This  choice 
linearizes  the  results  somewhat  and  eliminates 
a  considerable  amount  of  curvature  from  the 
plots. 

The  first  set  of  results  (fig  21) 
represents  the  data  for  a  nitrogen  pressure 


driver,  with  the  projectile  restraining  section 
between  gun  sections  1  and  2.  This  is  a  long 
barrel,  short  driver  case  with  a  relative  pin 
position  of  0.117.  Each  plus  sign  represents  an 
experimental  result.  Except  for  an  occasional 
"low"  point,  the  data  follow  Pidduck-Kent 
theory.  Referring  to  figure  17,  a  vertical  line 
drawn  at  a  relative  pin  position  of  0.1 17  would 
intersect  only  Pidduck-Kent  curves,  regardless 
of  initial  acceleration. 

Figure  22  shows  the  results  obtained 
with  a  nitrogen  driver  and  the  projectile 
restraining  section  between  gun  sections  3 
and  4.  The  relative  pin  position  is  0.372.  At  in¬ 
itial  accelerations  up  to  about  1000  g,  Pidduck- 
Kent  is  favored,  above  1000  g,  the  results  are 
closest  to  and  within  a  few  percentage  points 
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of  infinite  driver  theory.  Again  referring  to 
figure  17,  a  vertical  line  at  0.372  would  in¬ 
tersect  Pidduck-Kent  curves  up  to  800  g  and  in¬ 
finite  driver  curves  for  higher  initial  accelera¬ 
tions. 

The  third  set  of  data  (fig  23)  was 
based  on  a  configuration  that  was  geometrical¬ 
ly  the  same  as  the  last  set,  but  the  driver  gas 
used  was  helium.  In  this  set,  Lagrange  theory 
or  Pidduck-Kent  theory  applies  equally  to  the 
data  below  1000  g.  At  higher  acceleration,  the 
data  slightly  exceed  Pidduck-Kent  theory 
Reference  to  figure  16  at  a  relative  pin  position 
of  0.372  yields  only  Pidduck-Kent  intersections 
up  to  initial  accelerations  of  3000  g.  Figure  18 
shows  intersections  with  infinite  driver  only 
above  9000  g. 

The  combined  experimental  results 
are  presented  in  figure  24  for  comparison.  The 
data  of  each  set  (delta,  plus,  and  times  sym¬ 


bols)  are  sufficiently  separated  from  those  of 
the  other  sets  so  that  it  is  relatively  easy  to  pic¬ 
ture  the  trends  associated  with  each  Ex¬ 
perimental  results  at  about  200-g  initial  ac¬ 
celeration  show  equivalent  or  slightly  higher 
velocities  for  nitrogen  compared  with  those  for 
helium  at  the  same  projectile  restraining  sec¬ 
tion.  As  the  initial  acceleration  increases,  the 
superior  results  for  helium  are  apparent. 

The  experimental  velocities  attained 
with  nitrogen  at  a  relative  pin  position  of  0.372 
exceed  those  at  a  relative  pin  position  of  0.1 17 
(fig.  24).  However,  the  difference  in  velocities 
decreases  with  increasing  initial  acceleration. 
These  two  sets  of  data  should  eventually  cross 
when  the  0.1 17  gun  acceleration  is  sufficiently 
large  for  infinite  driver  theory  to  govern.  Then 
for  larger  initial  acceleration,  the  longer  barrel 
of  the  gun  with  a  relative  pin  position  of  0.117 
(0  883)  produces  higher  velocities  than  that 
with  0.372  (0.628). 
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Figure  21.  Data  and  theories,  nitrogen,  relative  pin  position  of  0.117. 


22 


VELOCITY  [FTfsJ  VELOCITY  IFIVsJ 


INITIAL  ACCELERATION  (g1/2) 

Figure  22.  Data  and  theories,  nitrogen,  relative  pin  position  of  0.372. 
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Figure  23  Data  and  theories,  helium,  relative  pin  position  of  0.372. 
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Figure  24  Combined  data,  nitrogen  and  helium,  relative  pin  positions  of  0.117  and  0.372. 


7.  DISCUSSION 

The  experimental  data  for  helium  show 
higher  velocities  than  for  nitrogen  at  the  same 
initial  accceleration  for  velocities  above  750 
ft/s  (225  m/s)  (fig.  24).  Helium  follows  Lagrange 
theory  or  Pidduck-Kent  theory  for  the  lower 
velocities  and  Pidduck-Kent  for  velocities 
above  1500  ft/s  (500  m/s)  (fig.  23).  From  the 
assumptions  associated  with  these  theories,  it 
follows  that  as  the  projectile  restraining  sec¬ 
tion  is  moved  away  from  the  muzzle  (decreas¬ 
ing  relative  pin  position),  these  theories  should 
continue  to  be  appropriate.  (It  is  not  clear  that 
these  theories  will  be  appropriate  if  the  pin  is 
shifted  toward  the  muzzle.  In  this  region,  the 
gun  approaches  an  infinite  driver  configura¬ 
tion,  and  infinite  driver  theory  may  govern  ) 

No  clear-cut  optimum  relative  pin  position 
has  resulted  from  the  experimental  or 


theoretical  investigation.  Choice  of  relative  pin 
position  becomes  a  matter  of  judgment,  which 
may  be  influenced  by  secondary  considera¬ 
tions  such  as  (1)  expected  range  of  use,  (2) 
minimization  of  the  volume  (mass)  of  gas  used, 
(3)  minimization  of  the  pressure  behind  the  pro¬ 
jectile  at  muzzle  exit,  and  (4)  ease  of  operation. 

Figure  25  presents  Lagrange  curves  and 
Pidduck-Kent  curves  for  the  first  three  possible 
projectile  restraining  section  locations  in  the 
gun  Analysis  of  Lagrange  theory  readily  yields 
the  initial  acceleration  required  at  one  projec¬ 
tile  restraining  section  compared  with  the  ac¬ 
celeration  at  a  second  projectile  restraining 
section  (acceleration  ratio)  to  produce  the 
same  muzzle  velocity.  (From  the  form  of  the 
Pidduck-Kent  equations,  it  appears  that  this  ac¬ 
celeration  ratio  varies  with  the  initial  accelera¬ 
tion.  Nevertheless,  from  the  Lagrange  curve 
and  Pidduck-Kent  curve  shapes  of  figure  25, 
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any  generalization  based  on  trends  applicable 
to  Lagrange  theory  is  generally  applicable  to 
Pidduck-Kent  theory  in  this  acceleration 
regime.) 

Consider  the  quantity  (mass)  of  gas  used 
per  test  as  a  secondary  criterion.  As  an  exam¬ 
ple,  assume  that  helium  is  used  and  that  the 
gun  is  operating  in  the  Lagrange  regime  or  the 
Pidduck-Kent  regime.  (If  the  acceleration  is 
sufficiently  high  to  operate  in  the  infinite  driver 
regime,  then  reducing  the  relative  pin  position 
increases  muzzle  velocity.)  The  gas  required 
for  a  shot  is  immediately  available  from  the  gun 
dimensions  and  required  initial  pressure  (ac¬ 
celeration)  (program  WTVSRPP). 

Figures  26  and  27  show  the  mass  of  gas 


(relative  to  the  projectile  mass)  required  versus 
projectile  restraining  section  location  (relative 
pin  position)  for  velocities  from  300  to  1000  ft/s 
(90  to  300  m/s).  A  significant  observation  is  that 
the  minimum  gas  mass  required  is  at  a  relative 
pin  position  of  zero.  (Although  as  the  relative 
pin  position  approaches  zero  the  required 
pressure  becomes  infinite,  the  volume  and  the 
mass  go  to  zero.)  The  quantity  of  gas  increases 
monotonically  with  increasing  relative  pin  posi¬ 
tion.  It  is  apparent  that  to  minimize  the  amount 
of  gas  used,  one  would  choose  the  smallest 
driver,  whether  it  be  helium  (fig.  26)  or  nitrogen 
(fig.  27).  This  choice  is  important  when  the  gun 
operates  with  helium  because  of  the  time  in¬ 
volved  in  purging  and  filling  the  driver;  that  is, 
the  smaller  the  driver  and  the  less  gas  re¬ 
quired,  the  faster  the  operation  (firing). 


Figure  25.  Lagrange  and  Pidduck-Kent  theories,  helium,  relative  pin  positions  of  0.117,  0.244,  and  0.372. 
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RATIO  OF  GAS  TO  PROJECTILE  MASS  RATIO  OF  GAS  TO  PROJECTILE  MASS 


RELATIVE  PIN  POSITION 

Figure  26.  Helium  (Lagrange  theory). 


RELATIVE  PIN  POSITION 

Figure  27.  Nitrogen  (Lagrange  theory). 
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Another  consideration  is  the  amount  of 
energy  stored  in  the  gas  when  the  projectile 
exits  from  the  gun.  This  energy  may  damage 
equipment  or  interfere  with  the  impact  ex¬ 
periments.  The  theory  presented  in  program 
NRG  predicts  the  energy  transferred  from  the 
driver  gas  to  the  projectile.  The  final  driver  gas 
energy  (enthalpy)  was  obtained  by  subtracting 
the  projectile  muzzle  energy  from  the  initial 
gas  energy  The  results  (fig.  28)  show  that  the 
smaller  the  driver  (small  relative  pin  position), 
the  greater  the  amount  of  energy  transferred  to 
the  projectile  and  therefore  the  smaller  the 
amount  of  energy  remaining  in  the  gas.  (Note: 
The  greater  the  degrees  of  freedom  of  the  gas, 
the  less  energy  transferred  to  the  projectile.) 

The  effects  of  projectile  restraining  sec¬ 
tion  placement  on  several  parameters  are 
presented  in  table  2.  The  five  projectile 
restraining  section  locations  are  nearest  the 


breech.  The  effects  are  arbitrarily  referenced 
to  the  present  pin  location  (between  sections  3 
and  4)  For  example,  the  initial  acceleration  re¬ 
quired  to  produce  a  particular  muzzle  velocity 
for  a  projectile  restraining  section  between 
sections  3  and  4  is  assigned  a  value  of  1.  The 
initial  acceleration  required  to  produce  the 
same  velocity  for  the  pin  between  sections  1 
and  2  would  be  2.02  times  as  great  (table  2, 
second  parameter). 

Note  that  the  fourth  pin  position  (relative 
pin  position  of  0.500)  is  close  to  the  optimum 
location  based  on  Lagrange  theory  (relative  pin 
position  of  0.465,  helium,  table  1).  Table  2  in¬ 
dicates  that  moving  the  pin  to  this  position 
reduces  the  initial  acceleration  by  only  3  per¬ 
cent,  but  requires  30  percent  more  gas.  If  a 
21 -percent  increase  in  initial  acceleration  is 
tolerable,  then  gas  usage  can  be  cut  21  per¬ 
cent  by  moving  the  projectile  restraining  sec- 
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tion  back  one  section.  Moving  the  projectile 
restraining  section  back  two  sections  to  its 
original  location  is  not  as  attractive. 

The  data  of  the  fourth  parameter  of  table  2 
show  that  there  would  be  a  24-percent  reduc¬ 
tion  in  gas  energy  at  breech  exit  by  moving  the 
projectile  restraining  section  back  one  section 

8.  CONCLUSIONS  AND  REMARKS 

A  theoretical  and  experimental  study  was 
made  to  determine  where  a  fixed  length,  con¬ 
stant  diameter  tube  should  be  divided  so  that 
when  one  portion  of  the  tube  is  used  as  a  gun 
barrel  and  the  remaining  portion  is  used  as  a 
driver  (chamber),  the  maximum  projectile  muz¬ 
zle  velocity  can  be  obtained  for  a  particular 
driver  gas  and  initial  acceleration.  Studies  of 
several  HDL  4-in.  gas  gun  configurations  in¬ 
dicated  that  no  one  configuration  is  best  for  all 
velocities  greater  than  500  ft/s  (150  m/s),  the 
range  of  interest.  (For  lower  velocities,  the  gun 
is  used  in  the  1-atm  configuration,  where  the 
entire  tube  length  is  used  for  the  barrel  and 
ambient  air  is  the  driver  This  configuration  is 
not  relevant  to  the  study  ) 

For  the  cases  considered  (short  driver  ver¬ 
sus  long  driver,  helium  versus  nitrogen,  and  in¬ 


itial  accelerations  from  200  to  5000  g),  the  ex¬ 
perimental  results  generally  follow  different 
theories  in  different  operating  regimes: 
Lagrange  theory  in  the  very  low  Mach  number 
regime,  infinite  driver  theory  in  the  very  high 
acceleration  regime,  and  Pidduck-Kent  theory 
in  other  regimes.  Although  the  various  results 
are  explicable  in  terms  of  the  postulated 
theories,  no  one  theory  accurately  predicts  the 
entire  range  of  experimental  results. 
Therefore,  it  is  difficult  to  infer  a  general  op¬ 
timum  configuration  because  theoretical 
predictions  around  the  optimum  are  in  most 
cases  conflicting  or  poorly  defined  or  both. 
Because  no  clear-cut  division  of  the  tube  can 
be  inferred  to  optimize  muzzle  velocity,  secon¬ 
dary  criteria  such  as  quantity  of  gas  required, 
energy  remaining  in  the  gas,  (statistical)  utiliza¬ 
tion  of  the  gun,  and  ease  of  operation  can  in¬ 
fluence  the  configuration  without  unduly  affec¬ 
ting  muzzle  velocity. 


In  conclusion,  significant  savings  in  gas 
usage  and  reductions  in  gas  energy  at  muzzle 
exit  can  be  realized  with  only  a  small  accelera¬ 
tion  penalty  by  reconfiguring  the  HDL  gun  to  a 
one  section  longer  barrel  and  one  section 
shorter  driver. 


TABLE  2.  EFFECTS  OF  PROJECTILE  RESTRAINING  SECTION  WITH  HELIUM  DRIVER 


Projectile  restraining  section  between  sections 
Parameter  _ 


1  and  2 

2  and  3 

3  and  4 

4  and  5 

5  and  6 

Relative  gun  position 

0  117 

0.245 

0.372 

0.500 

0  628 

Initial  acceleration  ratio  to  at¬ 
tain  particular  muzzle  velocity 

2  02 

1.21 

1  00 

0.97 

1.07 

Ratio  of  required  gas  masses 

0.63 

0.79 

1.00 

1.30 

1.81 

Ratio  of  final  driver  gas  energy 

0.46 

0  76 

1.00 

1.22 

1.42 

to  initial  gas  energy 


Note  Reference  is  projectile  restraining  section  location  between  sections  3  and  4;  relative  pin 
position  is  0  3  72. 
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NOMENCLATURE 

a  Initial  sound  speed  in  driver  reservoir 
A  Cross-sectional  area  of  gun  barrel 

K  2APL/M 

L  Total  gun  length  (barrel  +  driver,  X  +  Y) 
M  Projectile  mass 

N  Degrees  of  freedom  of  driver  gas:  N  =  3 
for  helium  and  N  =  5  for  nitrogen 
P  Initial  gas  pressure  in  driver  reservoir 
u  Projectile  Mach  number  based  on  initial 
sound  speed  (U/a) 

U  Projectile  velocity 

X  Length  of  driver  reservoir 

Y  Length  of  barrel 

Z  Relative  pin  position  in  gun  (X/L) 
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The  computer  programs  on  the  following  pages  were  written  in  VAX-1 1  BASIC,  a  compilable, 
native  mode  language.  Since  the  language  is  not  as  restrictive  as  standard  BASIC,  a  few  of  the 
major  differences  are  pointed  out  for  those  who  wish  to  adapt  the  programs  to  their  own  com¬ 
puters 

Comments. — In  addition  to  REM  statements,  the  exclamation  point  serves  as  a  comment 
delimiter. 

Example: 

10  THIS=THAT  !  SET  THE  VARIABLE  THIS  EQUAL  TO  THE  VARIABLE  THAT 

Continuations. — A  line  may  be  continued  on  the  following  line  when  the  last  character  on  the 
line  to  be  continued  is  an  ampersand. 

Example: 

10  THIS=  & 

THAT 

Multiline  statements. — Statements  to  be  executed  in  continuous  order  can  be  on  the  same 
line  if  they  are  separated  by  a  back  slash. 

Example. 

10  THIS=7  \THAT=8 

Line  numbers. — Line  numbers  are  not  required  if  the  statement  is  indented  at  least  one 
character  from  the  left  margin.  Lines  are  executed  in  order. 

Example: 

10  THIS=7 
THAT=8 

Variables. — Names  of  variables  are  allowed  to  be  up  to  32  characters  long  and  can  contain 
most  ASCII  characters. 

Example: 

10  DEG.OF.FREEDOM=2/(GAMMA-1)  \  AVG  A=SUM  A/ND 

Use  of  these  enhancements  produced  programs  that  are  easier  to  read  and  to  follow  than 
standard  BASIC.  The  programs  included  are  PDQ,  CONDIA,  CATCH,  CDPK,  THEORY,  WTVSRPP, 
and  NRG.  Other  programs  were  used  to  retrieve  the  data  from  the  record  files  and  to  plot  the 
data.  However,  these  are  not  general  and  would  be  of  little  use  to  the  reader. 
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10  REM' 


PDQ.BAS 

APkIL  1980  H  D  CURCHACK 


PIDDUCK-KENT  THEORY 
SEE  SEIGEL  HIGH  SPEED  GUNS  PP  180-181 


20  DIAM=4  !  IN.  ! 

GAS$(1)="HELIUM" 

SPEED( 1)=3300  ! FT/ SEC! 

DEG  OF  FREEDOM( 1 )=3 
GAMMA(T)=5/3 
THISSTEP=500 
TOTAL  L=94 


\  X  SEC  AREA=PI*DIAM'>2/4 

\  CSs$CZ)="nitrogen" 

\  SPEED(2)=1 145 

\  DEG  OF  FkEED0M(2)=5 

\  GAMHA(7)=7/5 

\  DEL  MU=1/THISSTEP 


30  FOR  GAS=1  TO  2  \  PRINT  GAS$(GAS) 

N0=DEG  OF  FkEEDOM(GAS)  \  CO=SPEED(GAS)  \  GAMA=GAMMA(GAS) 
GAS  STEp(T)=.0005  \  GAS  STEP(2)= .0025 


OPEN  "PK."+LEFT$(  GASS  (GAS)  ,3)  FOR  OUTPUT  AS  FILE  //GAS 
GOSUB  40  FOR  INIT  G=200  TO  3001  STEP  200 
CLOSE  II GAS 


OPEN  "PK.1  . "+LEFT$(GAS$(GAS)  ,3)  FOR  OUTPUT  AS  FILE  IlG AS 


GOSUB  40  FOR  INIT  G=3000  TO  15001  STEP  3000 
CLOSE  //GAS 
NEXT  GAS 
GOTO  90 


40 


!  PIDDUCK  KENT  SUBROUTINE 

I  - 


PRINT  "INITIAL  G  =";INIT_G  !  TERMINAL  OUTPUT  DURING  kUN 

!  A_BAR  IS  A  PARAMETER  USED  BY  SIEGEL  TO  DETERMINE  G/M  AND  VEL 

FOR  A_BAR=0  TO  1  STEP  GAS_STEP(GAS) 

!  INTEGRATION  ! 

INTGRAL=0 

FOR  INTVAL=0  TO  THISSTEP 
MU=INTVAL*DEL_MU 
F  OF  MU=( 1-A  BAR*MU*2)*(N0/2) 
n7TGEAL=  INTCKAL+F_OF_MU 
NEXT  INTVAL 

INTGRAL=( INTGRAL- . 5* ( 1+F  OF  MU) ) /THISSTEP 
!  INTEGRATION  COMPLETE  ! 


50  REM 


ALSO 


G_0VER_M=( NO+2) *A_BAR* ( 1-A_BAR) " ( -( NO+2) / 2) * INTGRAL 

-  CONVERT  G  OVER_M  TO  INIT  G  AND  DRIVER  L  - 
CO* 7=GAMA*P/RH0 

WEIGHT  OF  GAS  (G)=RHO*X  SEC  A*DRIVER  L 
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DRIVER_L=G/(RHO*X_SEC  A) 

=G*C0"2/ (GAMA*P*X  SEC  A) 

=G* CO" 2/ (GAMA* IN  IT  G*K) 

THEREFORE  “ 

60  DRIVER_L=G_OVER  M*C0"2/ (GAMA*INIT  G*32.2) 

FRAC  L=DRIVER  L7T0TAL  L  “ 

IF  FlAC_L>l  “  “ 

THEN  bO 

ELSE  VEL=2*C0/(GAMA-1)*SQR(A  BAR* (1 -FRAC  L"(GAMA-1) ) ) 
PRINT  IIG AS  USING  "//  .###F  ,  IHHHt . it"  ,THAC  L.VEL 
70  NEXT  A  BAR  - 

80  PRINT  ?GAS,"  1.  ,  0." 

RETURN 

90  END 
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10  REM' 


CONDIA 


H  D  CURCHACK 


APR  I960 


20  DIA=4  !  IN.  ! 


GAS$(1)=" HELIUM" 

SPEED( 1)=3300  ! FT/SEC! 

DEG  OF  FREEDOM( 1 )=3 
GAMHA(T)=5/3 
TOTAL  L=94 


\ 

\ 

\ 

\ 

\ 

\ 


X  SEC  AKEA=P I*DIA"2/ 4 

GT&if  (7)  NITROGEN” 

SPEED(2)=1145 

DEG  OF  FREEDOM(2)=5 

GAMMA(7)=7/5 

OU T$-“#.####  ,  ####.#" 


30  FOR  GAS=1  TO  2 

NO=DEG  OF_FREEDOM( GAS ) 

Nl=NO+T 

N01=N0/N1 

CO=SPEED(GAS) 

PRINT  GAS$(GAS) 

EXTEN$=LEFT$ ( GAS$ ( GAS ) ,3) 

F ILE$ ( 1 ) =" COND IA . "+EXTEN$ 

F ILE$ ( 2 ) =" BEG IHCD . "+EXTEN $ 

F ILE$ ( 3 ) =" CATCH . "+EXTE NS 

OPEN  FILE$(  I)  FOR  OUTPUT  AS  FILE  0GAS+I-1  FOR  1=1  To  2 
OPEN  FILE$(3)  FOR  INPUT  AS  FILE  #(G AS+2) 

GOSUB  40  FOR  INIT  G=0  TO  3000  STEP  200 
CLOSE  ilGAS  JGAS+I  JGAS+2 

FILE$  ( 1)="C0NDIA1  ,"-cEXTEN$ 

FILE$( 2)="BEGINCD1 . "+EXTEN$ 

F 1LE$ ( 3 ) =" CATCH 1 . "+EXTEN$ 

OPEN  FILE$(  I)  FOR  OUTPUT  AS  FILE  //GAS+I-1  FOR  1=1  TO  2 
OPEN  FILES (3)  FOR  INPUT  AS  FILE  il(G AS+2) 

GOSUB  40  FOR  INIT_G=0  TO  15000  STEP  3000 
CLOSE  #GAS  ,  //GAS+1 ,  //GAS+2 

NEXT  GAS 
GOTO  110 

40  !  INFINITE  RESERVOIR  SUBROUTINE 


PRINT  INIT  G 
IF  INIT  G<70 

THEN  PRINT  IlG AS,"  0.  ,  0." 

FOR  DRVR=0  TO  TOTAL  L  STEP  1 

G  X  L=32.2*lTOTAL  L-DRVR)  ! (FT/ SEC) SQUARE  ! 
GUSUB  60  “ 

PRINT  IIG AS  USING  OUTS  ,DRVR/TOTAL_L,UO*CO 
NEXT  DRVR 

50  INPUT  #( GAS+2) ,FRAC  L, 

G_X_L=  3  2 . 2*  ( 1  -FRACJE)  *  TOTAL_L 
GOSUB  60 

PRINT  HGAS+1)  USING  OUT$,FRAC  L,UO*CO 
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RETURN 

60  !  DETERMINE  VELOCITY  WHICH  IS  IMPLICIT  IN  TERMS  OF  DISTANCE. 

Y0=G  X  L*INIT  G/OT2 
Fl-.T  ~ 

U0=S0R(Y0) 

IF  U0>N0— 1 

THEN  U0=N0-I 
70  X9=Y0/ 10000 

IF  X9< 1 . 00000E-07 

THEN  X9=l .00000E-07 

!  LOOP 

X=N01  *  (1-  (1-U0/N01)  /  (1-U0/N0)^N1  ) 

X1=Y0-X 

S0=SGN(X1) 

IF  F2-S0 
THEN  F2=0 

Fl-Fl/2 

90  IF  ABS(X1)>X9 

THEN  U0=( 1+S0*(F1 ) )*U0 
F2=-S0 
GO  TO  80 

!  END  LOOP 

100  RETURN 

110  END 
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10  REM' 


CATCH 


APRIL,  1980 


CURCHACK 


Compute  point  at  which  first  reflected  wave  intercepts  the 
projectile  path. 

Elimination  of  the  velocity  at  which  the  first  wave  overtakes 
the  projectile  (Seigel  13-7  by  substitution  into  eq.  12-1) 
yields  the  distance  that  the  wave  overtakes  the  projectile. 


LET 

F« 1  -U*  (  GAMA- 1 )  /  (  2*  A  ) 

THEN 

F=( l+( GAMA+1 ) *X_D/ 2) * (-2* ( GAMA- 1) / (GAMA+1 ) )  (13-7) 

HENCE 

"  X  B=2/ (GAMA+1)* ( (2-( GAMA+1 )*F )/ (GAMA- 1) /F*( (GAMA+1) / (GAMA-1) )+l) 


WHERE  X  D  AND  X  B  ARE  THE  DIMENSIONLESS  DRIVER  AND  BARREL  LENGTHS 

\  X  SEC  AREA=PI*DIAM~2/4 


20  DIAM=4  !  IN.  ! 

TOTAL  L=94 
GAS$(T)  -"HELIUM" 

SPEED( 1)=3300  ! FT/SEC! 

DEG_0F  FREED0M( 1  )=3 
GAMA(lJ=5/3 


\  GAS$(2)="NITR0GEN" 

\  SPEED( 2)=1 145 

\  DEG  OF  FREEDOM(2)=5 

\  GAMA(  27=7/5 


30  FOR  GAS»1  TO  2 

N0=DEG  OF  FREEDOM(GAS) 

PRINT  "UAS7(GAS) 

OPEN  "CATCH. "+LEFT$(GAS$( GAS) ,3)  FOR  OUTPUT  AS  FILE  #GAS 
OPEN  "CATCH1  ,"+LEFT$(GAS$(GAS)  ,3)  FOR  OUTPUT  AS  FILE  ilG AS+2 
FI  GAM=(GAMA( GAS)+1) / 2 
F2  GAM=(GAMA(GAS)-l)/2 
C  1>Q=SPEED(GAS)“2 


FYL=GAS 

GOSUB  40  FOR  INIT_G=0  TO  3000  STEP  200 
CLOSE  #FYL 

FYL=GAS+2 

GOSUB  40  FOR  INIT_G=0  TO  15001  STEP  3000 
CLOSE  #FYL 
NEXT  GAS 
GOTO  90 


40  !  OVERTAKE  COMPUTATION 

1  - 

PRINT  "INITIAL  G  -";INIT_G  !  TERMINAL  OUTPUT  DURING  EXECUTION 

D  START=0 

D7STEP*TOTAL_L/8 

RESOLUTE*0 

TRAVEL  D-0 
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50  FOR  DRVR=D_START  TO  TOTAL  L  STEP  D  STEP 

XJD=32.2*DRVR*INl’f_G/C_Sl)  “ 

!  REFLECTED  WAVE  COMPUTATION 

F=( 1+F 1  GAM*X  D)“(-2*F2  GAM/FI  CAM) 

X_B=((1'=F1  GAff*F)/F“(Fl  GAM/F2  GAM)/F2  GAM+1)/F1  CAM 
IF  INIT  C<70  “  “  -  - 

THEN“TRAVEL  D-X  B*C  SQ/ (32. 2* INIT  G) 

60  CATCH  D=(TRAVEX  D+URVRj/ TOTAL  L  “ 

IF  CATCH  D<1  “  “ 

THEN  70 

ELSE  RESOLUTE=RESOLUTE+l 
IF  RESOLUTE>10 
THEN  SO 

ELSE  D_START=DRVR-D  STEP 
D_STEP=D  STEP/7 
GOTO  50 

70  NEXT  DkVR 

SO  FRAC  L=DRVR/ TOTAL  L 

PRINT  III  YL  USING  ,  #.####  t  ilHIHHt" ,  FRAC  L, CATCH  D,IN1T  G 


90  END 
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10  REM'' 


CDPK.BAS 


APk.  80  HD  CURGHACK 


20 


Read  previously  determined  Pidduck— Kent  velocity  from  the 
appropriate  file,  compute  the  infinite  driver  length  velocity 
at  the  same  relative  length,  then  choose  the  value  from  the 
theory  that  predicts  the  lower  velocity.  Once  the  latter 
theory  produces  the  lower  velocity ,  only  use  that  theory  since 


Pidduck-Kent  is  no  longer 

appropriate • 

DIAM-4 

\ 

X  SEC  AREA=PI*DIAM''2/4 

GAS$(1)=" HELIUM" 

\ 

GHSSCZ) -"NITROGEN" 

SPEED( 1)=3300  ! FT/SEC! 

\ 

SPEED(2)=1145 

DEG  OF  FREEDOM( 1)=3 

\ 

DEG  OF  FREEDOM(2) =5 

GAMMA(T)=5/3 

\ 

GAMHA(7)=7/5 

F$(l )-"#.###  ,  ####.#" 
TOTAL_L=94 

\ 

F$(2)«F$(1)+"  ,  *" 

30  FOR  GAS=1  TO  2 

NO-DEG  OF  FREEDOM (GAS)  \  PRINT  GAS$(GAS) 
fx$=le7t$Tgas$(gas) ,3) 

OPEN  "PK."+FX$  FOR  INPUT  AS  FILE  #(GAS+2) 

OPEN  "CDPK. "+FX$  FOR  OUTPUT  AS  FILE  II GAS 
OPEN  "CDPKINT . "+FX$  FOR  OUTPUT  AS  FILE  #(GAS+4) 

GOSUB  40  FOR  INIT  G-200  TO  3000  STEP  200 
CLOSE  #GAS,//(GAS+T)  ,#(GAS+4) 

OPEN  "PK1."+FX$  FOR  INPUT  AS  FILE  #(GAS+2) 

OPEN  "CDPK1 . "+FX$  FOR  OUTPUT  AS  FILE  #GAS 
OPEN  "CDPKINT1."+FXS  FOR  OUTPUT  AS  FILE  #(GAS+4) 

GOSUB  40  FOR  INIT  G=3000  TO  15001  STEP  3000 
CLOSE  #GAS,#(GAS+2M(GAS+4) 

NEXT  GAS 
GOTO  170 

40  PRINT  "INITIAL  G  =";INIT_G  !  TERMINAL  OP  DURING  RUN 

V  FLAG=0 

50  iHPUT  Z/GAS+2 ,FRAC  L,VEL 

DRVk=TOTAL  L*FRAC“L 

G  X  L-32. 2*\T0TAL~L-DRVR)  ! (FT/ SEC) SQUARE  ! 

GUSlJB  100 

NEW_VEL=UO*SPEED( GAS)  \  OUT_VEL=VEL 
IF  V  FLAG-1 

THEN  OUT  VEL-NEW  VEL 
GOTO  70 

60  IF  NEW  VEL<=VEL 

THE¥  OUT  VEL-NEW_VEL 

v  fTag= 1 

PHINT  #(GAS+4)  USING  F$(l),  FRAC_L ,OUT_VEL 
70  PRINT  IfGAS  USING  F$(  1+V_FLAG)  ,  FRAC_L  ,OUT_VEL 

80  IF  FRAC  L<1  THEN  50 
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90  RETURN 

100  !  SUBROUTINE  TO  DETERMINE  VELOCITY  GIVEN  A  LENGTH  FOR 

!  A  CONSTANT  DIAMETER  GUN.  EQUATION  IS 

!  IMPLICIT  FOR  VELOCITY  AND  EXPLICIT  FOR  LENGTH. 

Y0=G  X  L*1NIT  G/ SPEED(GAS) “2 

! - - - - - 

no  fi=  .  i 

1=0 

U0=SQR( YO) 

IF  U0>N0-1 

THEN  U0=N0-1 
120  X9=Y0/ 10000 

IF  X9<1 .00000E-07 

THEN  X9=l .00000E-07 

!  LOOP 

130  X=(N0/ (NO+1) )  *  (1-  (1-(N0+1)*U0/N0)  /  ( 1-U0/N0r  (NO+1)  ) 


!  SEE  FOR  EXAMPLE  EQ  B1 1 ,  HDL-TR-1330,  CURCHACK 

X1=Y0-X 

S0=SGN(X1) 

140  IF  F2=S0 

THEN  F2=0 

Fl=Fl/2 

150  IF  ABS(X1)>X9 

THEN  U0=(1+S0*(F1))*U0 
F2=-S0 
GO  TO  130 

160  RETURN  !  END  LOOP 

170  END 
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10  REM 


aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 


APR  80 


THEORY 


H  D  CURCHACK 


Range  of  Speeds  vs  Initial  g  for  Isentropic  Gun. 

Same  for  Constant  Dia  Gun.  Outputs  to  separate  files. 
Sane  for  Pidduck-Kent  Gun.  Outputs  to  separate  files. 


CONSTANTS 

20  DIA=4  !  IN.  \ 

SH0RT= 1  \  REG=2 

BARREL(SHORT)=59  \  BARREL(REG)=7 1 

DRIVER( SH0RT)=35  \  DRIVER(REG)=23 

RPP$( SH0RT)=" . SHO; 1"  \  RPP$(REG)=".MED; 1" 

G_X_L( I)=32.2*BARREL( I)  FOR  I=SH0RT  TO  LNG 
T0TAL_L=BARREL( SH0RT)+DRIVER( SHORT) 

GAS$(1)=" HELIUM"  \ 

SPEED( 1)=3300  ! FT/ SEC!  \ 

DEG  OF  FREED0M(1)=3  \ 

GaMMA(T)=5/3  \ 

F0RM2 $-"#.###  ,  #.###"  \ 

MAX_INIT_G( 1)=5500  \ 

30  FOR  GAS=1  TO  2 

S  T  ART_G= MAX_  I N IT_G( GAS ) 

CU=SPEED(GAS) 

GAMA=GAMhA( GAS ) 

NO=DEG_OF_FREEDOM( GAS ) 

M0RE$=LEFT$( GAS$ ( GAS) ,3) 

FOR  BARRELX=SHORT  TO  LNG 

MY_FILE$=" ISEN"+M0RE$+RPP$ ( BARRELS ) 

OPEN  MY  FILE$  FOR  OUTPUT  AS  FILE  ilGAS 
PRINT  'TILE=";MY_FILE$ 

FRAC  L*DRI VER( BARRE  L% ) / TOTAL  L 
FRACTL_FACT-NO*FRAC_L* ( 1-FKAC_L~ ( 2/ NO ) ) / 2 
FOR  INIT  G-0  TO  START  G  STEP  25 

YO-C  X  L(BARREL%)*INIT  G/ CO*  2 
K-2*XyTJ*T0TAL  L/BAkRELTBARRELX)) 
MACH_NO=SQR(  K*  FRAC_L_FACT ) 

PRINT  #GAS  USING  F0RM4$,INIT  G.MACH  NO*C0 
40  NEXT  INIT  G  -  “ 

CLOSE  #GA^ 

MY  FILE$-"COND"+MORE$+RPP$(BARREL%) 

OP'EN  MY  FILE$  FOR  OUTPUT  AS  FILE  VG AS 
PRINT  ,fFILE=";MY  FILE$ 

G  X  L(3)=»G  X  L(B7£RREL%)  !  (FT/ SEC) SQUARE  ! 

50  FOR  INIT  0=0  TO  START  G  STEP  25 

KO-INIT  G/X_SE‘C_AREA 
GOSUB  170 

PRINT  #GAS  USING  FORM4$ ,  INIT  G  ,  U0*C0 


X_SEC_AREA=PI*D1A''2/4 
\  LNG=3 

\  BARREL(LNG)=83 
\  DRIVER(LNG)=11 
\  RPP$(LNG)=" .LON; 1" 

!  (FT/ SEC) SQUARE 


GAS$(2) ="NITR0GEN" 
SPEED(2)=1145 
DEG  OF  FREEUOM(2)=5 
GAMHA(7)=7/5 

FORM4 $-"####.#  ,  ####.#" 
MAX  INIT  G(2)=2750 
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60  NEXT  INIT  G 

CLOSE  II GAS 
70  NEXT  BARRELS 

J  PIDDUCK-KENT  THEORY 

!  SEE  SEIGEL  HIGH  SPEED  GUNS  PP  180-181 

OPEN  "P IDD"+MORE$+RPP$( BRLS)  FOR  OUTPUT  AS  FILE  #BRL%  & 

FOR  BRLS=1  TO  3 

FLAG( 1)=0  \  FLAG(2)=0  \  FLAG(3)=0 
GAS_STEP( 1 )=. 001  \  GAS  STEP( 2)=.0025 
GOSUB  100  ~ 

CLOSE  III, 112, 113 
80  NEXT  GAS 
90  GOTO  220 


100  !  A_BAR  IS  A  PARAMETER  USED  BY  SEIGEL  TO  DETERMINE  G/M  AND  VEL  t 
THISSTEP=500  \  DEL  MU=1/THISSTEP 
FOR  A_BAR=0  TO  1  STEP  GAS  STEP(GAS) 

FOR  SUMMY=0  TO  THISSTEP-1- 
MU=SUMMY*DEL  MU 
NUMU=MU+DEL_MU 

DEL_AREA=  (  1  -A_BAR*MU ~  2 ) ~ ( NO/ 2 ) +( 1 -A  B  AR*  NUMU'' 2 ) ~ ( NO/ 2 ) 
INTGRAL=INTGRAL+DEL  AREA  ~ 

NEXT  SUMMY  ~ 


INTGRAL=INTGRAL/2/THISSTEP 

G  OVER  M=(N0+2)*A  BAR*(1-A  BAR)~(-(N0+2)/2)*INTGRAL 
FUR  BRE%=SHORT  TO- LNG  — 

FRAG  L=DRI VER( BRLS ) /TOTAL  L 
IF  F7aG(BRL%)=1  ~ 

THEN  130  !  DON'T  COMPUTE  IF  INIT  G  TOO  LARGE 

INIT  G(BRLS)  =  G  OVER  M*C0“'2/ (GAMA*DRlVER  L*32.2) 
VEL=7*C0/(GAMA-iy*SQRTA  BAR* (1 -FRAG  L"(GAMA-1))) 

IF  INIT_G(BRLS)<=START  G  ~ 

THEN  PRINT  #BRLS  USTNG  FORM4$,INIT  G( BRLS). VEL 
ELSE  FLAG(BRLS)=1  ~ 

NEXT  BRLS 
INTGRAL-0 

IF  INIT_G(1)>STAR1_G 
THEN  RETURN 
150  NEXT  A_BAR 
160  RETURN 


110 


120 


130 

140 


!  COMPUTATION  SUBROUTINE 

170  YO-G  X  L(3)*R0*X  SEC  AREA/C0~2 

F1=*.T 

1-0 

UO-SQR(YO) 

IF  UO>NO— 1 

THEN  U0=N0-1 
180  X9-YO/ 10000 

IF  X9< 1 . 00000E-07 

THEN  X9-1.00000E-07 
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!  LOOP 

190  X=(N0/ (NOt-1)  )  *  (1-  (1-(N0+1)*U0/N0)  /  ( 1-U0/N0r  (NO+1)  ) 

1=1+1 
X1=Y0-X 
S0=SGN(X1) 

IF  F2=S0 
THEN  F2=0 

F 1 =F 1 / 2 

200  IF  ABS(X1)>X9 

THEN  U0=( 1+SO* (FI ) )*U0 
F2=-S0 
GO  TO  190 

!  END  LOOP 

210  RETURN 
220  END 
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10  REM 

ix 

6. 


AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 


WTVSRPP.BAS 


6c 


6c 


The  isentropic  expansion  theory  permits  ready  estimation  of 
the  quantity  of  gas  (proportional  to  G.OVLk.M)  required 
to  produce  a  given  velocity  as  a  function  of  the  relative 
pin  position  (FRAC.L).  Note  that  G.OVER.M  is  a  monotonic 
function  of  FRAC.L. 

-  CONVERT  G.OVER.M  TO  INIT.G  AND  DRIVER. L  - 
C0"2=GAMA*P/RHO 

ALSO 

WEIGHT  OF  GAS  (G)=RH0*X. SEC . A*DR1VER. L 
DRIVER . L=G/ ( RH0*X. SEC . A) 

=G*C0"2/ (GAMA*P*X. SEC. A) 

-G*C0"2/(GAMA*INIT.G*M) 

THEREFORE 

DRIVER. L=G. OVER. M*C0"2/ (GAMA* INIT.G*32. 2) 

15  DIM  FR( 200) ,GOM(B ,200) 

20  GAS$( 1)=" HELIUM"  \  GAS$( 2) “"NITROGEN" 

GAMMA(l)=5/3  \  GAMMA(2)=7/5 

SPEED(1)=3300  ! FT/ SEC!  \  SPLLD(2)=1145 

fi  §="#.////#  it. mm  it  .it  it  it  it  if  it  it.itititm  it.  it  it  it  it  it  it"  +  & 

"  it.  it  it  it  it  it  if  it.  it  it  it  it  it  it  it.  it  it  it  it  it  it  it.  it  it  it  it  it  it" 

30  FOR  GAS-1  TO  2 

PRINT  GAS$(GAS); 

GAMA-GAMMA(GAS) 

G AMI -GAMA- 1 

GAM2  -GAMA* GAM  1/2  \  PRINT  "FACTOR  GAM2  =";GAM2 

MAX .  L= .  8 

K-0 

FOR  VELOCITY-300  TO  1000  STEP  100 
K-K+l 
KK=0 

MACH.NO— VELOCITY/ SPEED( GAS) 

F ACT 1 -MACH. NO" 2  *  GAM2  \  PRINT  VELOCITY, FACT1 
FOR  FRAC.L-0  TO  MAX.L  STEP  .005 

G.0VER.M-FACT1  /  (  1-FRAC . L"GA111) 

IF  G.OVER.M  <=1 
THEN  KK-KK+1 

. G0M( K , KK) -G.OVER.M 
FR(KK) -FRAC.L 
ELSE  GOTO  37 

35  NEXT  FRAC.L 
37  NEXT  VELOCITY 

OPEN  "WTVSRPP0."+LEFTS(GAS$(GAS),3)  FOR  OUTPUT  AS  FILE  ItG AS 

MARGIN  #GAS , 100 

PRINT  #GAS  USING  Fl$ ,FR( J) ,GOM(8 , J) ,G0M( 7 , J) ,GOM(6 , J) ,  6. 

GOM(5 ,J) ,GOM(4 , J) ,G0M(3 , J) ,G0M( 2 , J) ,G0M( 1 ,  J)  FOR  J=1  TO  KK 
CLOSE  #GAS 
NEXT  GAS 
40  END 
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10 


REM . 


NRG. BAS 


AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 


Computation  of  fraction  of  gas  energy  converted  to  projectile 
energy  based  on  Lagrange  theory.  Function  only  of  Z  (or  RPP) 
and  U  (degrees  of  freedom). 


20 

30 

35 

40 

50 

60 


Using  gas  energy  per  unit  mass,  e 
and 

and  total  gas  energy, 
and 
Then 
or 

• 

Also 

and  from  eql  of  text 
Therefore,  Eb 


mass,  e 

= 

N  R  T  /  2 

P 

= 

RHO  R  T 

Lgas 

= 

e  RHO  Volume 

volume 

= 

A  x. 

hgas 

(N  R  T  /  2)  X  RHO  A  x 

Egas 

N  P  A  x  /  2. 

2 

Ebird 

= 

MU  /  2 

2/N 

Ebird 

= 

Egas  (1  -  Z  ) 

2/N 

i  /  Egas 

= 

1  -  Z 

DEF  FNA(Z ,N)=( 1-Z~ ( 2/N) ) 

FORM $-"#.####  #.#####  il.iHHHHI" 

OPEN  "ENERGY.DAT"  FOR  OUTPUT  AS  FILE  #1 

PRINT  #1  USING  F0RM$,Z,FNA(Z,3) ,FNA(Z,5)  FOR  Z=0  TO  .049  STEP  .005 
PRINT  #1  USING  FORM$,Z,FNA(Z,3),FNA(Z,5)  FOR  Z=.05  TO  1  STEP  .025 
END 
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